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Abstract 

The algebraic method provides useful techniques to identify models in 
designs and to understand aliasing of polynomial models. The present 
note surveys the topic of Grobner bases in experimental design and then 
describes the notion of confounding and the algebraic fan of a design. The 
ideas are illustrated with a variety of design examples ranging from Latin 
squares to screening designs. 

Contents 

1 Ideals and varieties: introducing algebra 

2 Grobner bases 

2.1 Term orderings 

2.2 Matrix based term orderings 

2.3 Monomial ideals and Hilbert series 

2.4 Grobner bases 

2.5 Software tools 

3 Experimental design 

4 Understanding aliasing 

5 Indicator functions and orthogonality 

6 Fans, state polytopes and linear aberration 

6.1 The algebraic fan of a design 

6.2 State polytope and linear aberration .... 

7 Other topics and references 

1 Ideals and varieties: introducing alge- 
bra 

We are familiar with the use of polynomials throughout statistics. For 
example, much of this handbook is concerned with design for polynomial 
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regression. Thus we have polynomial terms in factors x\,x%,...: 

xi, x'i, XtXi, . . . 

and a second order polynomial response surface in two factors is 

f{xi,X2) = 0oo + O w xx + 6qix 2 + 2 ox\ + O11X1X2 + o2 x 2 2 . (1) 

The first, but very important, algebraic point is that polynomials are 
made up of linear combinations of monomials. Consider a set of k factors 
xi, . . . ,Xk and non-negative integers a — (ai, . . . , ecu)', a monomial is 

a a 1 ac a k- 
x — X X ' ' X 

Note that when we use the term polynomial we shall typically mean a 
polynomial in one or more variables. 

A monomial x a can be represented by its exponent vector a and we 
can list the monomials in a model either directly or by listing a set of 
exponents. We shall often use the notation {x a ,a £ M}, for some set 
of exponents, M. This chapter is largely concerned with the interaction 
between the choice of a design and the list M. We know from classical 
factorial design that only some models are estimable for a given design and 
so any such theory must be intimately related to the problem of aliasing 
and we shall cover this is section [4] 

The set of all polynomials over a base field is a ring, so that rings are 
the basis of the theory. Thus, given a base field K we obtain the ring of 
polynomials, R = K[xi, . . . ,xt] over K, which are linear combinations of 
monomials with coefficients in the base field. Our "a" notation allows us 
to write this compactly as 

/O) = eaX ™> 

where, as above, M is a finite set of distinct exponents and clearly f(x) £ 
R. For example, the set M for the polynomial model in Equation {TJ is 
{(0,0), (1,0), (0,1), (2,0), (1,1), (0,2)}. 

Given that we have launched into algebra we need to introduce the 
first two essentials: ideals and varieties. In what follows we present only 
the basic ideas of the theory, pointing the reader to [11] or [41] for further 
details. 

For a ring R we have special subsets called ideals. 

Definition 1 A subset I C R is an ideal if for any f,g £ I we have 
f + g £ I and for any f € I and g £ R we have fg £ J. 

The ideal generated by a finite set of polynomials . . . , f m } is the set 
of all polynomial combinations: 

{h, ■ ■ ■ , fm) = H h f m g m ■ gi, ■ ■ ■ ,g m e R} 

To have some immediate intuition consider a single point x. The set 
of all polynomials / such f(x) = is an ideal: since for any polynomial 
g(x) if f(x) = we have g(x)f(x) — 0. In the next section this will be 
extended to sets of points, namely designs. The Hilbert basis theorem 
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says that any (polynomial) ideal / is finitely generated, i.e. for any ideal 
we can find a finite collection fi,...,f m G R such that / = {/i, . . . , f m ). 

We are familiar with linear varieties expressed by setting some linear 
polynomial function equal to zero. Thus a straight line can be written as 
the collection of points (2:1,2:2) in two dimensions such that axi + bx2 + c — 
0, for constants a, b, c. An algebraic variety is the extension of this concept 
to simultaneous solutions of a set of polynomial equations. 

Definition 2 Let /1, . . . , f m G K[xi, . . . , x k ] be a set of polynomials. The 
associated affine variety is the solution ( also called zero set ) of a set of 
simultaneous equations they define: 

V(fi,.- .,/„) = {(01,.. - ,o fc ) G K k : fi(ai,.. .,a d ) = 0, i = l,...,m} 

Every afhne variety has an associated ideal which we write I(V). It is the 
set of all polynomial which are zero on the variety: 

I(V) = {f€K[ ] : f(ai, . . . , cifc) = 0, for all (ai, . . . , a k ) G V}. 

What appears to be a straightforward relationship between ideals and 
varieties is actually very subtle. If we start with polynomials /1 , . . . , f m 
and construct the corresponding variety V and form the ideal I{V), is it 
true that I(V) = . . . , /m)? We can always claim that . . . , f m ) C 
/(V), but the converse may not be true and refer to [TT] for a detailed 
discussion. Fortunately, for a design, the variety is collection of isolated 
single points, the equivalence holds and we may move freely between ideals 
and designs. 

2 Grobner bases 

Perhaps the most important construction in abstract algebra is that of a 
quotient. Give two polynomials f,g€ K[xi, . . . , Xk] and an ideal I define 
the equivalence class / ~/ g if and only if / — g 6 /. The members of the 
quotient K[xi, . . .x„]/I are the equivalence classes. Since /1 ~j fi and 
gi ~i g% imply fi+gi ~i /l+ffa and fagi ~j f%g%, then K[xi, . . .x„]/I is 
also a ring. Finding K[x\, . . . x„]/I in a particular case requires a division 
algorithm. Finding a quotient computationally needs a division algorithm. 

2.1 Term orderings 

Let us recall division of polynomials in one dimension. If we divide 1 + 
3a; + 2a; 2 + x 3 by 2 + x we would obtain the tableau: 

a; 2 + 3 

a; + 2 ) x 3 + 2x 2 + 3x + 1 
-a; 3 - 2x 2 

3a; + 1 
—3a; — 6 

giving: x A + 2x 2 + 3x + 1 = (a; 2 + 3) (a; + 2) — 5. We give this example to 
remind ourselves that at each stage we need to use the leading term. To 
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obtain leading term we need an ordering. In one dimension the ordering 
is 1 -< x -< x 2 -< • • • , That is, we order by degree and division is unique. 
This is generalised to a special total ordering on monomials {a; 1 *}- 

Definition 3 A monomial term ordering, -<, is a total ordering of mono- 
mials such that 1 -< x a for nil a > 0, a / and, for all j > 0, x a -< x 13 
implies x a+ ~' -< a; /3+7 . 

We shall use the term monomial ordering for short. There is a number 
of standard monomials orderings. 

1. Lexicographic ordering, Lex. x a -<Lex x 13 when (i) /3 — a > and the 
leftmost entry of ft — a is positive. 

2. Graded lexicographic ordering, DegLex: x a -<DegLex x 13 if (i) the 
degree of a is less than that of /3, jaj < \/3\ and (ii) a <Lex P 

3. Reverse lexicographic ordering, DegRevLex: x a -<De 3 revLex x^ if (i) 
\a\ < /3\ and (ii) a -<Lex P, where the overline means: reverse the 
entries. 

Graded orderings are orderings in which the first comparison between 
monomials is determined by their total degree. For example, under a 
graded order, x\ y Xj for any indeterminates x%, Xj in the ring K[x\, . . . , a;*] 
The degree lexicographic and degree reverse lexicographic term orders 
above fall in this class. Contrary to graded orderings, for a lexical or- 
dering in which Xi >- Xj then Xi >- x™ for m = 1,2,... thus making all 
powers of Xj lower than x%. 

2.2 Matrix based term orderings 

Monomial term orderings can be defined using products with matrices and 
element- wise comparisons. If the exponents of monomials x a ,x® are con- 
sidered as row vectors, we say that x a -<m x 13 if Ma T < Mf} T ', where M 
is a non-singular matrix and the inequality is tested element-wise starting 
from the first element. The matrix M above satisfies certain conditions 
which are stated in the following theorem [41| . 

Theorem 4 Let M be a full rank matrix of size k x k such that the first 
non-zero entry in each column is positive. Then M defines a term ordering 
in the following sense: 

1. For every vector a £ Z> with a/(0,...,0) then (0, . . . , 0) < Ma T 
and 

2. For any pair of vectors a, fi, 7 G Z> such that Ma T < Mf5 T then 
M(a + 7) T < M(P + j) T . 

The identity matrix of size k corresponds to the lexical term ordering. 
Note that the relation between ordering matrices and term orderings is 
not a one to one. A matrix M' defining the same ordering as M can 
be obtained by multiplying each row of M by a positive constant so for 
instance the matrix with diagonal 1,2,..., k and zeroes elsewhere also 
defines a lexical term ordering. Usually only integer entries are used for 
computations although the theory does not preclude using for instance, 
matrices with rational or real entries |11] , 
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An important case of ordering matrices is that of matrices for graded 
orderings. Any full rank matrix M in which all elements of the first row 
are a positive constant defines a graded ordering. The degree lexicographic 
term ordering is built with a matrix M with all entries one in its first row 
and the remaining rows are the top k — 1 rows of an identity matrix. The 
CbCbA command Use T : : =Q [x , y , z] , DegLex ; creates the same ring and 
ordering when the matrix and ring are defined with the commands 

M:=Mat([ [1,1,1] , [1 , 0, 0] , [0, 1 ,0] ] ) ; 

Use T: :=Q[x,y,z] , Ord(M) ; 

The querie xy~2>x~2z; yields output FALSE which means that xy 2 -< 
x 2 z under the graded lexicographic order in which x y y y z. 

The standard ordering in the software system CbCbA is the degree 
reverse lexicographic (DegRevLex), which is implicit in the following ring 
definition 

Use T: :=Q [x,y,z] ; 

xy~2>x~2z; 

The output of the querie is TRUE and this is interpreted as xy 2 y x 2 z 
under a degree reverse lexicographic term ordering in which x y y y 
z. Note the reversal of the ordering between the two monomials for the 
previous graded order. 

A more specialized and efficient instance of matrix orderings is pro- 
duced by a using a single row matrix, in which case we say "ordering 
vector". An ordering vector defines only a partial but not a total order- 
ing over R. For example the vector w — (1, 1, 1) naturally produces the 
ordering xy 2 y w xz because 

(M,l)(l,2,0) T = 3>2 = (l,l,l)(l,0,lf, 

yet it cannot distinguish between monomials of the same degree such as 
xy 2 and z 3 . However, Grobner basis are computed over finite sets of 
monomials rather than over all monomials with exponents in Z> . This 
last fact together with a careful selection of the ordering vector are at the 
core of the efficient Universal Grobner bases algorithms [Tl 131]. 

2.3 Monomial ideals and Hilbert series 

Now that we have a total ordering any finite set of monomials has a 
leading term. In particular, since a polynomial, /, is based on a finite set 
of monomials it has a unique leading term. We write it LT^(f), or, if -< 
is assumed, just LT(f). 

A monomial ideal is an ideal generated by monomials. Monomial ideals 
play a critical part in computational methods for polynomials. 

Definition 5 A monomial ideal I is an ideal for which a collection of 
monomials f\ , . . . , f m such that any g £ I can be expressed as a sum 

m 

g = ^2gi(x)fi(x). 
1=1 

Multiplication of monomials is just achieved by adding exponents: 

a B a+R 

Ob Ob 0Z 
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and a + f3 is in the positive (shorthand for non-negative) "orthant" with 
corner at a. The set of all monomials in a monomial ideal is the union of 
all positive orthants whose corners are given by the exponent vectors of 
the generating monomial fi, ■ ■ ■ , f m - 

For a given monomial ideal, a complete degree by degree description 
of the monomials inside the ideal or, equivalently, those outside the ideal 
is given by the Hilbert function and series. Here we only give the basic 
idea, referring the reader to references [TT] and [12] for a full description. 

Definition 6 Let I be a monomial ideal in R. 

1. For all non-negative degrees j, the Hilbert function HFi(j) is the 
number of monomials not in I of total degree s. 

2. the Hilbert series of I is the formal series HSi(s) = ~}2JL Q s 3 HF(j). 

The Hilbert series is the generating function of the Hilbert function. Both 
count monomials which are not in the monomial ideal 7. In what follows, 
unless it is required, we omit the subindex referring to the monomial ideal. 

Example 1 Consider the ideal 7 = {x 3 ,xy 2 ,y 4 } C k[x,y]. The monomi- 
als which do not belong to 7 are 1, x, x 2 , y, yx, yx 2 , y 2 , y 3 , so the Hilbert 
function equals 1, 2, 3, 2 for j = 0, 1, 2, 3 and zero for all j > 4. The Hilbert 
series is thus HS(s) = 1 + 2s + 3s 2 + 2s 3 . 

Note that monomials in the first orthant are counted with the formal 
series X^jlo ( J+ ^~ 1 ) s "' w h ere k is the number of indeterminates. Then, by 
substracting the Hilbert series HS(s) from the last expression we have a 
generating function to count monomials inside 7. Using two dimensions 
as in Example [1] the formal series for the firtst orthant is 



oo / . \ oo 



which counts the monomials in the first quadrant. The generating function 
for the number of monomials in 7 for each degree is found by substraction: 



(1 + 2s + 3s 2 + 2s 3 



2s 3 + s 4 



(l-s) 2 ' ' ' ; (l-s) 2 ' 

The alternating signs in the polynomial in the numerator are related to 
inclusion-exclusion rules and although it may seem a simple calculation 
above, in general determining a simple form of the numerator in this last 
computation is not a simple task. 

Example 2 Consider the monomial ideal in k[xi, . . . ,xj] generated by 
monomials x 2 , . . . , x 2 and by all pairs XiXj, 1 < i < j > < 7. This ideal has 
Hilbert Function with values 1 and 7 for j = 0, 1 and zero for j > 2 so its 
Hilbert Series is HS(s) = 1 + 7s, i.e. one monomial of degree zero and 
seven monomials of degree one outside the ideal. The monomials outside 
this ideal are 1, xi, X2, £3, X4, x$, x§, x-j, which later will be understood as 
a model for the Plackett-Burman design in Examples [6l and 1111 See also 
first row of Table 0] 
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2.4 Grobner bases 



Dickson's Lemma states that, even if we define a monomial ideal with 
an infinite set of fi, we can find a finite set h\,...h n such that / = 
{hi, . . . ,hk). But there are, in general, many ways to express a ideal I as 
being generated from a basis I — {/i, . . . ,/ m ). 

Definition 7 Given an ideal I a set {g±, . . . g m } is called a Grobner basis 
if: 

(LT{ 9l ),...,LT{g m )) = <LT(/)>, 

where (LT(I)) is the ideal generated by all the monomials in I. 

We sometimes refer to (LT(I)) as the leading term ideal. 

Lemma 8 Any ideal I has a Grobner basis and any Grobner basis in the 
ideal is a basis of the ideal. 

Monomial orderings are critical in establishing that for any given 
monomial ordering, -<, any ideal / has a unique "reduced" Grobner ba- 
sis. Given a monomial ordering, -<, and an ideal expressed in terms of 
the G-basis, I = (gi,...,g m ) with respect to that monomial ordering 
any polynomial / has a unique remainder, r(x) with respect the quotient 
operation K [xi, . . . , Xk]/I- That is 

m 

/ = X>(x)Si(aO + r(a:) (2) 

i = l 

We call the remainder r(x) the normal form of / with respect to / and 
write r(x) — NF(f). Or, to stress the fact that it may depend on -<, we 
write NF(f,~<). 

The division of a polynomial in Equation @ is the generalization of 
simple polynomial division such as that of the example shown in Page O 
where the result was si(x) — x 2 + 3 with remainder r = —5. In other 
words the normal form of 1 + 3x + 2x 2 + x s with respect to the ideal 
generated by gi(x) — 2 + x is —5. 

Here are some formal definitions. 

Definition 9 Given a monomial ordering -<, a polynomial f — ~^2 aeL S a x a 
is a normal form with respect to -< if x a £ (LT(f)) for all a £ L. 

Lemma 10 Given an ideal I and a monomial ordering -<, for every f £ 
K[xi, . . . ,Xk] there is a unique normal form NF(f) such that f—NF(f) £ 
I. 

We now need to relate (i) the Grobner basis, (ii) a division algorithm 
and (iii) the nature of the normal form. We have partly covered this but 
let us collect the results together. 

1. There are algorithms, which given an ideal, /, a monomial ordering 
-< and a polynomial / deliver the remainder r, in Lemma (8), by 
successively dividing by the G-basis terms gi, i = l,...,m. The 
best known is the Buchberger algorithm. 

2. Suppose the remainder r(x) — NF(f) = YlaeL 6aX a , then {x a , a £ 
L} is precisely the set of monomials not divisible by any of the leading 
terms of the G-basis of / : LT(gi), i = 1, . . . , m. 
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3. The remainder r(x) = NF(f) does not depend on which order the 
G-basis terms Qi(x) are used in the division algorithm. 

4. The (maximal) set {x a ,a £ L}, which can appear in a remainder 
r(x) is a basis of the quotient ring, considered as a vector space of 
functions over k\x\, . . . , x^]/I. The terms are linearly independent 
over I: 

a 

implies 9 a — for all a € L. 

2.5 Software tools 

All the operations denned above are available on modern computer algebra 
software. Here is a brief list, the full list is very extensive and extends 
to nearly all areas of computer algebra, sometimes called computational 
algebraic geometry: |9]CoCbA, [22], macaulay2, [28 gfan, 23 Singular. 
A rough list of capabilities relevant to this chapter is as follows. 

1. Construction of monomial orderings; the standard ones are usually 
named and immediately available 

2. Ideal operations such as unions, intersections, elimination 

3. Buchberger algorithm and modern improvements, quotienting, Nor- 
mal Forms 

4. Special algorithms for ideals of point. We shall use these extensively 
in our examples 

5. Grobner fan. See section [6] 

3 Experimental design 

We have indicated already that for applications to design we should think 
of design as lists of points, 

D = {x W ,...x (n) }, 

in R k . As algebraic varieties they have associated ideal 

1(D) = {/ : f(x) = 0, x e D} 

The use of polynomials to define design is clearly not new. For example 
a 2 k full factorial designs give by {±1, . . . , ±1} is expressed the solution 
of the simultaneous equations: 

{xi - 1 = 0, i = l,...,k}. 

To obtain fractions we impose additional equations: e.g. xi . . . Xk = 1. 

We now give what can loosely be described as the algebraic method in 
the title of this Chapter. We do this in a step-by-step approach. 

1. Choose a design D 

2. Select a monomial term ordering, -< 
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3. Compute Grobner basis for 1(D) for given monomial ordering, -<. 

4. The quotient ring 

K[x 1 ,...,x k ]/I{D) 

of the ring of polynomials K[xi, . . . , Xk] in xi, ■ ■ ■ , Xk forms is a vec- 
tor space spanned by a special set of monomials: x a ,a £ L. These 
are all the monomials not divisible by the leading terms of the G- 
basis and \L\ = \D\. 

5. The set of multi-indices L has the "order ideal" property: a 6 L 
implies f3 £ L for any < f3 < a. For example, if x\x2 in the model 

SO is 1, Xl, X2,X\X2- 

6. Any function y(x) on D has a unique polynomial interpolator given 
by 

f( x ) = 9aX °' 

such that y(x) — f(x), x £ D. 

7. The cardinality of the design and the quotient basis is the same: 
\L\ = \D\. 

8. The X-matrix is n x n, has full rank n and has rows indexed by the 
design points and columns indexed by the basis: 

X = {x a }igD,a£i 

The implications of the method are considerable. But at its most basic 
it says that we can always find a saturated polynomial f(x) interpolating 
data over an arbitrary design D. 

The shape of the model index set L arising from the order ideal prop- 
erty is important. It is exactly the shape which, in the literature has 
been called variously: "staircase models", "hierarchical models", "well- 
formulated models" , or "marginality condition" , see |37l 140] . It can be 
be seen easily from the fact that the multi-index terms given by L are 
the complement in the non-negative integer orthant of those given by the 
monomials in the monomial ideal of leading terms: the complement of a 
union of orthants has the staircase property. 

We now give a number of examples. 

Example 3 Screening designs. A class of designs for main effect esti- 
mation while simultaneously avoiding biases caused by the presence of 
second order effects and avoid confounding of any pair of second order 
effects was recently proposed |29] , The authors produced designs of size 
n — 2k + 1 for different dimensions ranging from k — 4 up to k = 30, and 
their construction is based on folding a certain small fraction of size k of 
a 3 fc design with levels —1, 0, 1 and then adding the origin. Naturally that 
after folding and adding the origin, the screening design still remains a 
special fraction of 3 fe design. Here we consider the designs for k = 7 and 
k = 10 in Tabled] 

For k = 7, the design is obtained by first folding the points 0+-+-+-, 
-0+-++-, +-0++++, + — 0+ — , — ++0 — , -+-++0+, +++++-0 and then adding 
the origin to total 15 points. Under the usual degree reverse lexicographic 



9 

















Xl 


X2 


X3 


X'4 


X5 


X6 


X7 


Xs 


Xg 


alio 



















1 


1 


1 


1 


1 


1 


1 


1 


1 



















— 1 


— 1 


— 1 


— 1 


— 1 


— 1 


— 1 


— 1 


— 1 


Xl 


X2 


X3 


X4 


xs 


xe 


Xl 


1 





— 1 


— 1 


— 1 


— 1 


1 


1 


1 


1 





1 


— 1 


1 


—1 


1 


— 1 


— 1 





1 


1 


1 


1 


—1 


—1 


— 1 


— 1 





— 1 


1 


— 1 


1 


— 1 


1 


1 


— 1 





— 1 


1 


I 


—1 


— 1 


1 


1 


— 1 





1 


— 1 


1 


1 


— 1 


— 1 


1 





1 


— 1 


— 1 


1 


1 


— 1 


— 1 


1 





— 1 


1 


— 1 


— 1 


1 


1 


— 1 


— 1 





1 


1 


1 


1 


— 1 


— 1 


1 


— 1 





1 


1 


1 


1 


— 1 


1 


1 





— 1 


— 1 


—1 


—1 


1 


1 


— 1 


1 





— 1 


—1 


—1 


—1 


1 


— 1 


1 


1 





— 1 


—1 


1 


— 1 


1 


1 


—1 


— 1 





1 


—1 


—1 


—1 


1 


— 1 


—1 





1 


1 


—1 


1 


—1 


-1 


1 


1 





-1 


1 


1 


1 


-1 


1 


1 


-1 





1 


-1 


1 


-1 


-1 


-1 


1 


1 





-1 


-1 


-1 


1 


-1 


-1 


1 





-1 


1 


-1 


1 


1 


1 


-1 


-1 





1 


1 


1 


1 


-1 


1 


-1 


1 





-1 


-1 


1 


-1 


1 


-1 


1 


1 





1 


-1 


-1 


1 


-1 


1 


-1 





1 


1 


-1 


1 


-1 


1 


-1 


-1 





-1 


1 


1 


-1 


1 


1 


-1 


-1 





1 


-1 


1 


1 


1 


1 


1 


-1 





-1 


-1 


1 


-1 


-1 


1 


1 





-1 


1 


-1 


-1 


-1 


-1 


-1 


1 





1 


1 


1 


-1 


-1 


1 


-1 


1 





-1 























-1 


-1 


-1 


1 


1 


-1 


1 


-1 





1 
















1 


1 


1 


-1 


1 


-1 


1 


-1 


-1 



















-1 


-1 


-1 


1 


-1 


1 


-1 


1 


1 


















































Table 1: Two screening designs [29 . 



ordering in CbCbA, we identify the model with terms: 1, xi,Xa,xa, xa, x$,X6, 
X7,X6,X7,X2X?,X3X7,X4X7,X5X7 and X6X7. We note that use of a graded 
order allows for the inclusion of all terms of degree one before the addition 
of terms of second degree, and the total degree of this model (addition 
of all exponents) is 21. If a degree lexicographic order is used, the model 
remains with the same total degree but it interchanges one interaction for a 
quadratic term: 1, x\, £2, 23, X4, xs,xe, xr, x\, x\, X7, X5X6, X4X7, X5X7, x%xt. 

Lexical term orderings work in rather the opposite manner than graded 
orderings. For a lexical ordering, then selection of terms is concentrated in 
including all terms with Xj. As this cannot go further than 1, 057, asf, then 
term selection allocates all possible terms including xq until interaction 
with x§X7 appears and term X6X7 can be allocated with xr, returning to xq 
again when required, eventually including terms with X5. The model has 
terms 1, xr, x r ,Xe, x 6 x 7 , x^xj, x%, x%X7, x\x^, x 5 , x 5 x 7 , x 5 x 6 ,x$x 6 X7, x\,x\xr 
and its total degree is 31. 

Varying term orders over all possible orderings is in general a complex 
and expensive task. In Section [6] we discuss and comment on the whole 
set of models identified by this seven factor design, when considering all 
possible term orders. 

For k = 10 the points to be folded are 0+++++++++, +0 ++++, 

+-0-++ — ++, + — 0++++ — , +-++0 — +-+, +-++-0+-+-, ++-+-+0 — +, ++-++ — 0+-, 
+++ — +-+0- and +++-+-+ — 0. The standard term ordering in CbCbA was 
used to identify with this design a model of total degree 30 with constant, 
all ten linear terms x\, . . . , asio, two quadratic terms x%, x\ and eight dou- 
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ble interactions between xio and one of X2, ■ ■ ■ ,Xg. A lexical term ordering 
produces a model of much higher total degree (44) which contains monomi- 
als in four factors only (x§, Xs, xg and a^io): 1, £10, iio, £9, xgxio, Xgx^, Xq, 
XqXio, :rgx? , ^8, £8£io, £8^9, £82:92:10, x%, SgXio, xe, 2:62:10, x^xg, xaxgxio, 
xqXs and X6XsXio 

Example 4 Response surface design, non-standard. Here we take a 16- 
point design which is a 5 2 factorial with all internal points (a 3 2 design) 
removed: 

(2,0) (2,1) (2,2) (1,2) 

(0,2) (-1,2) (-2,2) (-2,1) 
(-2,0) (-2,-1) (-2,-2) (-1,-2) 
(0,-2) (1,-2) (2,-2) (2,-1) 

The design ideal is generated by the following polynomials x\ — 5x% + 
Ax2,x\ — 5x\-\-Axi,x\x\ — Ax\ — 4:X% + l&. It can be shown that for any term 
ordering, the above polynomials form a reduced Grobner basis and thus 
the design identifies a single model with terms 1, X2, x%, x%, x%, x%, X1X2, 
x\x\tX\x\,x\x\ix\,x\x2,x\, x\x2, xf, x\x2- In Example 1151 a standard 
response surface design of the central composite type is presented. 

Example 5 Regular fraction. Let us take the resolution III in six vari- 
ables (all main effects estimated independently interaction). In classi- 
cal notation this has defining contrasts: {ABCD, CDEF} . Instead of 
A, . . . , F we use indeterminates x\, . . . , x§ and selecting one of the four 
blocks expressed we have the ideal 

(x\ — \,X% — 1, X% — \,x\ — 1, x\ — 1,X%— 1, XxX2X3Xa — 1, X S XaX 5 X & — 1), 

and setting all polynomials above equal to zero (simultaneously) gives the 
design. The design ideal is created in the following CbCbA code as the 
sum of the ideal defining the full factorial design and the ideal defining 
the desired fraction. 

Use T: :=Q[x[l. .6]] ; 
I:=Ideal([A~2-l|A In IndetsQ]) 

+Ideal (x [1] *x [2] *x [3] *x [4] -1 , x [3] *x [4] *x [5] *x [6] -1) ; 

The CbCbA command QuotientBasis (I) ; gives the quotient basis 

[1, x[6], x[5], x[5]x[6], x[4], x[4]x[6], x[3], x[3]x[6], x[2], 
x [2] x [6] , x [2] x [4] , x [2] x [4] x [6] , x [1] , x [1] x [6] , x [1] x [4] , 
x[l]x[4]x[6]] 

If the confounding relation is desired for a given monomial, this is com- 
puted using the normal form. For example NF(x [2] *x [3] *x [6] , I) ; with 
output x[l]x[4]x[6] shows that over the design, the term X2X3X6 is 
aliased with X1X4X6, equivalently X2XzXq — xix^x^ £ I{D) and thus 
both terms appear in the same row of the aliasing Table The alias- 
ing table is read row-wise e.g. the first row implies that over the design 
1 = x\X2X3Xi — xsx^xsXfj = X1X2X5X6. Note that the first column of 
Table [2] corresponds to the quotient basis computed before, and that the 
row containing the monomial 1 has the generators of the defining contrast. 

For regular fractions like this case, the effect of different term orderings 
in the model means selecting a (possibly) different representative per each 
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1 


X1X2X3X4 


X3X4X^X & 


X1X2X5XQ 


Xi 


X2X3X4 


X1X3X4X5X5 


X2X5X6 


X2 


XIX3X4 


X 2 X 3 X4X 5 X 6 


XlX 5 X 6 


X3 


XIX2X4 


X4X5XQ 


XiX2X 3 X 5 X 6 


X4 


X1X2X3 


X3X5X6 


XIX2X4X5XQ 


X 5 


XiX2X 3 X 4 X 5 


X3X4XG 


XIX2XQ 


x 6 


XiX2X 3 X4X 6 


X3X4X5 


XiX 2 X 5 


X\X4 


X2X3 


X1X3X5X6 


X2X4X5XQ 


X1X6 


X2X3X4XQ 


X1X3X4X5 


X2X5 


X2X4 


X1X3 


X2X3X5X6 




x 2 x 6 


XiX 3 X4Xe 


X2X3X4XS 


XlX 5 


X 3 X 6 


XIX2X4XQ 


X4X5 


XlX 2 X 3 X 5 


X4XQ 


X1X2X3XQ 


X3X5 


XIX2X4X5 


X 5 X 6 


XxX 2 X3X4Xf,XQ 


X3X4 


X1X2 


XIX4X6 


X2X3X6 


X1X3X5 


x 2 x 4 x 5 


X2X4X6 


X1X3XQ 


X2X3X5 


XIX4X5 



Table 2: Aliasing table for Example [5] 



row of the aliasing table. If for instance the ring is defined instead with a 
lexical ordering with the command Use T: :=Q[x[l. .6]] , Lex; then the 
model terms enter in a lexical fashion. Ten terms 1, X2, X4, £5, xg, X2X4, 
X2X6, X4X6, X5X6, X2X4X6 of the model coincide with the model identi- 
fied above and six terms xi, X3, X1X4, X1X6, X3X6, X1X4X6 are replaced by 
£2X5X6, X4X5X6, X2X4X5X6, X2X5, X4X5, X2X4X5- In each case, the replace- 
ment monomial is in the same row. 

Example 6 Plackett-Burman, PB(8). Consider the Plackett-Burman 
design [35] with 8 points in k = 7 dimensions generated by circular 
shifts of the generator +++-+ — together with the point +++++++. With 
the standard ordering in CbCbA, we retrieve the usual first order model: 
1, xi,X2,X3, X4, Xs, Xe, X7. If a lexical term ordering in which xi >- • • • >- X7 
is used, the model retrieved is a "slack" model in only four factors with 
terms 1, xa, X5, xq, X7, X5X6, X5X7, XQX7. 

Example 7 Latin Square. It is a straightforward exercise to code up 
combinatorial a designs using indicator variables. Let us take as an ex- 
ample the 4x4 Graeco-Latin square derived via the standard Galois field 
method. The square is 

Aa BP C7 D5 
Bj AS Da CP 
CS Dj Ap Ba 
DP Ca BS A-y 

Coding up the design with 0— 1 indicators: Xij, i, j = 1, . . . ,4, where i 
indexes the factors rows, columns, and Latin, Greek letters and j the fac- 
tor "levels". The design points in this coding are shown in Table [3] Using 
the a graded lexicographic term ordering in C0C0A, the model identified 
for the design is 
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[1, u[4], u[3], u[2], t[4], t[3], t[2], c [4] , c [3] , c [2] , r [4] , 
r[3], r[2], t[4]u[4], t[4]u[3], t[4]u[2]] 

where the factors labelled u,t identify treatments (Latin and Greek let- 
ters) and the factors labelled c,r identify rows and columns of the design. 
Note the neat decomposition of model terms that coincides with the stan- 
dard analysis of variance for this orthogonal design: 



Source 


d.o.f. 


Mean 


1 


u (treatment factor 1) 


3 


t (treatment factor 2) 


3 


r (row factor) 


3 


c (column factor) 


3 


interaction tu (error) 


3 


Total 


16 



The interaction between treatment factors (three terms involving tu above) 
is often allocated to the error. 
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Table 3: Design points for Graeco-latin design of Example [71 

Example 8 Balanced Incomplete Block Design, BIBD. Consider the bal- 
anced incomplete block design with n = 12 runs and t — 6 treatments 
1 1 , . . . , ta arranged in b — 6 blocks of size two [U , tj] for the following pairs 

(1,2), (2,3), (3, 4), (4,5), (5,6), (1,6). 

Using the standard term ordering in CoCoA gives the following model: 

[1, t[6], t[5], t[4], t[3], t[2], b[6], b[6]t[6], b[5], b[4], 
b[3] , b[2]] 
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A similar decomposition to that of Example would allocate the inter- 
action bate to the residual error with only one degree of freedom. Under 
a lexical ordering we retrieve the same model as above. This result is not 
extremely surprising given the highly restricted range of monomial terms 
for the model for this design. Thus the biggest influence in selection of 
model terms is given by the ordering of the indeterminates, also known 
as initial ordering, see [41] . 

Example 9 Latin Hypercube Sample. Latin hypercubes [31] are widely 
used schemes in the design and analysis of computer experiments. The 
design region is often the hypercube [0, l] fe and designs of interest are 
often those that efficiently cover the design region. Latin hypercubes 
have at least two clear advantages: univariate projections of the design 
are uniform and they are simple to generate. 

The design Li with points (0, 0, 0), (1/5, 1, 4/5), (2/5, 3/5, 2/5), 
(3/5,4/5,1/5), (4/5,1/5,1) and (1,2/5,3/5) is an example of randomly 
generated latin hypercube in k = 3 dimensions and n = 6 runs. Under 
the standard term ordering in CbCbA, the design L\ identifies the model 
1, xi, X2, X3, 3:22:3, x%. Experimentally, some latin hypercubes have been 
found to identify certain types of models which are of minimal degree 
called "corner cut models", see [35] also [1]. The design L\ belongs to 
such class, and will be discussed further in Section [6] 

A second example of latin hypercube is L2 with points (0,0,4/5), 

(1/5,1/5,2/5), (2/5,2/5,1), (3/5, 3/5, 0), (4/5, 4/5, 3/5) and (1,1,1/5). 

Under the same ordering as above, L2 identifies the model 1,X2, X3, X2X3, 

2 ^1 
x 3 ,x 3 . 

4 Understanding aliasing 

The algebraic method is not only a way of obtaining candidate models 
but it does, we claim, deliver considerable understanding of the notion of 
aliasing. Aliasing is close to the idea of equivalence used above to define 
the quotient operation. 

Let 1(D) be the design ideal and for two polynomials /, g define 

f(x) ~ c g(x) 

to mean f(x) = g(x), x £ D. This is equivalent to 

f(x)-g(x)eI(D). 

Again equivalently we have, with respect to a particular monomial order- 
ing -<, 

NF(f) = NF(g) 

We call this algebraic aliasing. 

However, this is not quite the same as the statistical idea of aliasing. It 
would be enough that / = eg over the design for some non-zero constant c. 
That is both / and g should not both be in the same regression model. We 
first need a notation to refer to values of a polynomial f(x) on the design 
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expressed as a vector we write this as supp D (/(a;)). Then f(x) g(x) 
is equivalent to 

supp D (/(z)) = supp D (g(x)) 

Definition 11 Collections of polynomials F and G are said to be statis- 
tically aliased if 

span{supp(f), f 6 F} = span{supp(g), g G G} (3) 

and let us write this as 

F^ D G 

Given that f(x) — NF(/), x G D, we can rewrite[3]as 

span{supp(iVf (/)), / £ F} = span{supp(iVf g G G}. 

This means that any aliasing statement is equivalent to one for the normal 
forms. For / 6 F, let 

QgL 

where L is as defined above and depends on the design D and the mono- 
mial ordering, -<. Let Of be the vector of a j and define 9 a , g , similarly. 
Then since the matrix X in non-singular, by construction, we have 

-F ~d G <4> span{6>/,/ G F} = span{6» 9 ,g G F} 

Thus, statistical aliasing can be thought of in two stages: (i) first re- 
duce to expressing each polynomial in F and G to its normal form using 
the algebra then (ii) compare the coefficient subspaces. In the regular 
factorial fraction case the normal form of a monomial is itself a mono- 
mial, which makes the interpretation easier, but in the general case it is 
a polynomial. 

We can often we can find the alias classes by inspection, once we have 
the normal form. Consider Example 2 and the monomials {x\x\, x\x 2 , x\x%, 
x\x%}. Then, using CbCbA the normal forms are, respectively, 

4x 2 + 4y 2 - 16, 16x 4 + 16y 4 - 256, 
320x 4 + 320y 4 - 256x 2 - 256y 2 - 4096, 
5376a; 4 + 5376y 4 - 5120a; 2 - 5120i/ 2 - 65536 

We see by inspection that 

d 2 2 4 4i ri 6 6 8 81 

{l,X 1 X 2 ,X 1 X 2 } ~D {l,X 1 X 2 ,X 1 X 2 \- 

The equivalence continues to all {1, x\ k x 2 2 , x\ k+1 xl k+l }. 

To retain the link to classical notation we might say that the collection 
{I, A 2 B 2 , A i B i } is aliased with the collection {I , A 6 B 6 , A 8 B 8 } and we 
might write 

{I, A 2 B 2 , A 4 B 4 } « {I,A 6 B 6 , A 8 B 8 } 

This arises because A 2 B 2 ^A 2 + B 2 -4I and A 4 B 4 w A 4 + B 4 - 16, and 
both the reduced forms are estimable. 
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In this example odd terms also pair up. The normal forms of {xfxt, %iX<2 



a?2, ajfa?!} are respectively 



4a; 3 y + 4xy 3 - lGxy, 80x 3 y + 80xy 3 - 384xy 
134Ax 3 y + 1344xy 3 - 6656xy, 21760a; 3 j/ + 21760a;i/ 3 - 108544xy 

So that {1, x\x\, x\x\} « {1, x\x\, x\x^\, and so on, and in classical no- 
tation: {I, A 3 B 3 , A 5 B 3 } « {/, A 7 B 7 , A 9 B 9 }. 



5 Indicator functions and orthogonality 

At times it is convenient to see the design D as a subset of a full factorial 
design A/". This is most usual when we start with some basic design, such 
as a full factorial, and consider a fraction. We saw such a fraction in the 
last subsection. In this case an algebraic description of the fraction is via 
an indicator function: Fd, rather than a G-basis. The design ideal of D is 
unique, what changes are the generating equations we choose to describe 
it. These encode different information on D. 

An indicator function is a single additional function which we add to 
the generators of the ideal of the full factorial design to form the ideal of 
D. We can write the last example as 

7(D) = (x\ - l,x% - l,x% - l,x 1 x 2 x 3 + 1). 

The first three terms form the G-basis of the full factorial {(±1, ±1, ±1)}. 

From the equation f(xi,X2,xs) = X1X2X3 + 1 = we can deduce the 
indicator functions of D in M as g(xi,X2,xs) = = |(— 2:12:22:3 + 1). 
This takes the value 1 on the design and on J\f \ D. Then, on D: 

X1X2X3 + 1 = <S> g(x 1 ,x 2 ,x 3 ) = 1 

More generally let A/" be the basic, starting design which is not neces- 
sarily a full factorial design, and let D C N be a fraction. Fix a monomial 
order and, via the I{N), construct a vector space basis for interpolation 
over A/". Then the indicator function of D interpolates the 0, 1 values as 
required: 

t \ T l, x £ D 

9(x> ~\ 0, x e M\D 

In the example above, there is only one basis for interpolation, being 
N a full factorial design: {2:" 1 2:2 2 2:3 3 : on £ {0,1} for i = 1,2,3} and 
the indicator function involves only the terms for a = (0, 0, 0) and a = 
(1,1,1). 

The coefficients of the indicator functions expressed over the interpola- 
tion basis embed information on the 'geometric/combinatoric' properties 
of the fraction. We exemplify this in the binary case where JV is the 2 d 
with coding { — 1, 1} [20]. For factors with mixed levels a coding with 
complex numbers is needed [43] . 

Two (square- free) monomials x a ,x^ are said to be are said to be or- 
thogonal over D C Af if the corresponding columns in the A-matrix are 
orthogonal: 

= E x<x+B = °- 

xeD' x£D' 
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We can express this in terms of the indicator function over {-1, l} d and 
write 

J2 * a+p g{*) = o 

because g(x) = over N \ D and g(x) = 1 over D. In the example 
above we want to check that the two-way factors are not orthogonal to 
the one-way factor. Indeed 

^ 2 — X\X 2 X 3 sr-^ 2XiX 2 X 3 — 1 ^ 

2 y asia^iCa = o = ^ ^i^a^ - 4 = -4 ^ 

aJ6.IV" aSAf ~ x£Af 

because x\ — 1 over A/". It is no coincidence that the coefficient of X1X2X3 
in the indicator function is not zero. Out of the zero coefficients of g one 
can deduce the orthogonal (monomial) functions over D. 

A very practical advantage of the indicator function is that we can 
take union and intersections of design rather easily by using Boolean type 
operations over D: 

3r>inD 2 = gD 1 gD 2 , gD l uD 2 = gn x + go 2 — gD 1 gD 2 - 

Again the zero coefficients of the normal form of gD 1 nD 2 an d 9d 1 ud 2 over 
the interpolation monomial basis of J\f are informative of the geometry of 
the intersection and union design. 

6 Fans, state poly topes and linear aber- 
ration 

The computations of Grobner basis and model identification with Grobner 
basis described in Sections [5] and [3] depend upon the term ordering se- 
lected. Setting a fixed term order allows the experimenter to put pref- 
erence over terms which will be identified by the model, for instance a 
graded ordering will include as many terms of order one as possible fac- 
tors before adding terms of second degree in the model. In other instances, 
the experimenter might be interested in exploring the range of all models 
identifiable by the design using algebraic techniques. For example this 
would allow assessment of design properties like estimation capacity [8] 
or the minimal linear aberration of the design [4] and its general case 
of non-linear aberration [3]. Fan computations have been applied among 
others, to industrial experiments [27] and systems biology [15] . 

Some figures of this Section were generated with gfan and computa- 
tions were performed with CoCoA and gfan [91 128]. 

6.1 The algebraic fan of a design 

Given a design ideal 1(D) and ranging over all possible term orderings, 
we have a collection of reduced Grobner bases for 1(D). A crucial fact is 
that despite the infinite number of different term orderings (excluding the 
trivial case of one dimension), this collection of bases has always a finite 
number of distinct elements [3S]. Associated to this collection of Grobner 
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Figure 1: Algebraic fan of design L\ of Example [9] 



bases there is a collection of polyhedral cones, called the Grobner fan, and 
we term the algebraic fan of the design to the collection of different bases 
for the quotient ring R/I(D). Note that the algebraic fan is effectively, a 
collection of saturated models. 

For some relatively simple designs, such as factorial designs, the al- 
gebraic fan has only a single model. The general class of designs with a 
single model is called echelon designs [3T]. However, at present, compu- 
tation of the algebraic fan of a design remains an expensive computation. 
Reverse search techniques are at the core of state-of-the-art software gf an 
[28] . However, other approaches remain under investigation, such as the 
polynomial-time approach based on partial orderings, operations with ma- 
trices and zonotopes [T] [3T] . The well known link between Grobner basis 
calculations and linear algebra operations for zero dimensional ideals (i.e. 
design ideals) allows these methodologies to be efficient [131 130] 

Example 10 (Continuation of Example O The collection of all mod- 
els identifiable by the design L\ (algebraic fan of Li) was computed. 
Design L\ identifies 27 different models which can be classified in only 
six types of models, up to permutations of variables: l,x\,x\,x\,x\,x\ 
(3 models); 1, xi, x\, x\ , x\ , X2 (6 models); 1, iei, x\ , a:? , X2, x\Xi (6 mod- 
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els); 1,xi,X2,Xi,xiX2,X2 (3 models); 1, x\, X2, xz, x\, x\ (3 models) and 
\,x\,X2,X3,x\,xiX2 (6 models). 

We say that this fan has a complete combinatorial structure, mean- 
ing that each class of models is closed under permutations of indetermi- 
nates, e.g. if the model 1, the fan, so are the models 

1, X2, X2, X2, x\ and 1, X3,x%, xl,x^,x^, obtained by permuting inde- 
terminates. 

The algebraic fan of L\ is depicted in Figure [T] where each model is 
represented as a staircase diagram, with indeterminates x\,X2,x$ along 
axes and one small box for each monomial term. The models are presented 
by classes following the order described above (row- wise from top left). 
For instance, the first diagram shows the model 1, xs, x%, x%, x 3 , x 3 , the 
second is 1, xi,x\,x\, xf, x\ and so on. 

Now we turn our attention to the other latin hypercube L-z- From 
the design coordinates we note that this design has complete confounding 
between x\ and X2 and we should expect a much more limited collection 
of models. Indeed this design identifies only 11 models which are depicted 
in Figure [2] Only one of the models contains terms with x\ (first from 
left in second row); while the rest of the models have monomials in X2 
and £3. The models can be classified in three classes, only one of which 
is closed under permutation of indeterminates (shown in the left column 
in Figure [2]). 




Figure 2: Algebraic fan of design Li of Example [9] 

Example 11 (Continuation of Example[6]) In total there are 610 different 
hierachical models identifiable by the Plackett-Burman design. Those 
models belong to 4 different classes, only two of which are generated by 
all permutations of factors. Note that as the design has only two levels in 
each factor, the models identified by this design are all multilinear. The 
lowest total degree of models is 7, and the largest total degree is 10. See 
Table 0] for details and examples for each class, where the sign * refers to 
a class which is closed under permutation of indeterminates. The Hilbert 
Series HS(s) has been included to describe model terms degree by degree 
in each class. 
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Simplicial 
complex 


Degree 


HS(s) 


Class size 


Example 
Vertex/Model 


. • • 


7 


1 - 


h7s 


1* 


(1,1,1,1,1,1,1) 

l,Xi,X2,X3,X4,X5,X6,X7 




8 


1 - 


h 6s + s 2 


105* 


(1,1,1,2,1,2,0) 

1, Xi, x 2 , x 3 , x 4 , x 5 , x 6 , x 4 x 6 


n 


9 


1 - 


h 5s + 2s 2 


420 


(0,1,2,2,3,0,1) 

1, X 2 , X 4 , X 3 , X 5 , X7, X 4 X 5 , X3X5 




10 


1- 


h 4s + 3s 2 


84 


(1,3,0,3,0,3,0) 

1, Xi, X2, X 4 , X6, X 2 X 4 , X2X6, X 4 X6 



Table 4: Summary of the algebraic fan of the Plackett-Burman design. 



Example 12 (Continuation of example [5| The fan of the regular 2 6 ~ 2 
fraction with generators {ABCD, CDEF} is of relatively modest size: 132 
models which belong to six equivalence classes whose size range from 12 to 
24. Models range from total degree 26 to 32 and none of the equivalence 
classes is closed under permutation of indeterminates, which is not entirely 
surprising given the regularity of the design. Despite this apparent fan 
simplicity, this six classes share only three different total degrees and 
Hilbert functions. For instance, three different model classes share the 
same total degree 26 while other two different model classes have total 
degree 28. Table [5] shows a summary of the fan computations for this 
design, and Figure [3] shows simplicial representation of models in each 
class (vertices refer to single factors, edges to two factor interactions and 
so on). 

U' n£ Q & ^ & 

Figure 3: Depiction of simplicial models for fan classes I- VI (left to right), design 
2 6 - 2 . 

Example 13 (Continuation of Example[3| The algebraic fan of the screen- 
ing design for seven factors k — 7 and n — 15 runs is a complicated and 
large object which nevertheless exhibits in some instances combinatorial 
symmetry. The design identifies 18368 staircase models which can be clas- 
sified in 25 equivalence classes. The class sizes range from 7 to 2520, while 
total degree of models range from 21 to 31. Six equivalence classes are 
closed under permutation of indeterminates, and this includes the classes 
of models identified by degree lexicographic (420 models) and by degree 
reverse lexicographic (210 models); examples of models for each order- 
ing were computed in Example [3] Other 2 equivalence classes are almost 
closed, each can be paired with a other small equivalence class. 
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Class J° tal HS(a) 
degree 



Class 



Example 
Vertex/Model 



II 26 



III 



1 + 6s + 7s 2 
+2s 3 



(5,5,6,2,1,7) 
24 1,xi,x 2 ,x 3 ,X4,x 5 ,X6,xix 2 ,xix 3 ,xiX6, 

X2X3, x 2 x e ,x 3 x 6 ,x 4: x e ,x 1 x 3 x 6 , x 2 x 3 x 6 

(7,1,4,4,3,7) 
24 l,xi,x 2 ,x 3 ,X4,x 5 ,x 6 ,xix 3 ,x 1 x 4: ,xix 5 , 

xix 6 , x 3 x 6 , X4,x 6 , x 5 x e , xix 3 x 6 , xiX4,x 6 

(4,4,6,2,2,8) 

24 l,X 1 ,X2,X 3 ,X 4 ,X a ,X e ,XiX 3 ,X 1 X 6 ,X 2 X 3 , 

x 2 x 6 , x 3 x 6 , X4x e , x a x 6 , XiX 3 x 6 , x 2 x 3 x 6 



IV 



V 



9fi 1 + 5s + 7s 2 
28 +3s 3 



(6,6,2,6,0,8) 
24 1, cci, X2, X3, X4, xq,xix 2 ,xiX4, xixe, x 2 X4, 

X 2 X 6 , X 3 X 6 , X4X 6 ,XiX 2 X(i, XiX4Xq, X 2 X 4 X 6 

(4,8,0,8,4,4) 

24 1, X 1 ,X 2 ,X4, X 5 ,X 6 ,XiX 2 ,X 1 X4, X 2 X4, X 2 X 5 , 

X 2 Xq, X4X5, X4Xq,X\X 2 X4, X 2 X4X 5 , X 2 X4Xq 



VI 



32 



1 + 4s + 6s 2 
+4s 3 + s 4 



(8,0,8,8,0,8) 
12 1, xi, x 3 , X4, x e ,xix 3 ,xiX4, xix e , x 3 x 4 , x 3 x e , 
X4Xe,xix 3 X4, xia;4a;6, a^i^3^6, x 3 X4Xe,xix 3 X4Xe 



Table 5: Summary of the algebraic fan of regular fraction 2 6 



6.2 State polytope and linear aberration 

The state polytope of 1(D) is a geometric object which is associated with 
the Grobner fan of 1(D) [2] [35]. The state polytope is constructed as 
the convex hull of state vectors, and each state vector is built from a 
model in the algebraic fan by simply adding the exponents of the model. 
Aside from a constant, indeed each state vector is the centroid of the 
staircase diagram represented by the model and thus the state polytope 
is the convex hull of all those centroids. 

The state polytope of 1(D) encodes information by variables about 
the total degree of each model in the fan of design D. A simple argument 
of linear programming shows that models in the algebraic fan are those 
that minimise a simple linear cost function on the weighted degree of the 
model. This is the idea of linear aberration defined in [4]. This concept 
has been generalised to nonlinear cost functions [3l 113] . 

Example 14 For the latin hypercube design L\ , the state polytope of its 
design ideal is built with state vectors for each of the 27 models ennumer- 
ated in Example 1101 For instance, the model 1 1 has state 

vector 

(0, 0, 0) + (1, 0, 0) + (2, 0, 0) + . . . + (5, 0, 0) = (15, 0, 0) 

and as the other two models in this class are created by permutations of 
variables, the same action is performed on the state vectors so for this 
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class we have three vectors: (15, 0, 0), (0, 15, 0) and (0,0,15). A similar 
construction and arguments are used for each model in the fan of L\ and 
we have 6 vectors with permutations of each of (10,1,0), (7,2,0) and 
(4, 2, 1) ; three permutations for each of (4, 4, 0) and (6, 1, 1). 

There is a special type of polynomial models which are of minimal 
weighted degree. These models are termed corner cut staircases [38], 
as their exponents can be separated by their complement by a single 
hyperplane. The properties of corner cut staircases and their cardinality 
have been studied in literature [101 148] . 

A design that identifies all corner cut models is termed a generic design, 
and automatically a generic design is of minimal linear aberration [1] . The 
collection of models identified by design L\ (of Examples [9] [10] and I14[) 
is the set of all corner cut staircases for k = 3, n — 6 and thus L\ is a 
generic design. State polytopes associated with corner cuts and generic 
designs were described in [36] . 




Figure 4: Grobner fan for designs L\ (left) and L2 (right) of Example [9] 

In addition to information about degrees of models in the fan, the 
state polytope also encodes information to compute Grobner bases. To 
each vertex of the state polytope, a normal cone is associated [49]. The 
collection of all those cones is precisely the Grobner fan of 1(D), in the 
sense that the interior of each full dimensional cone contains ordering 
vectors necessary to compute the Grobner basis (and identify the model) 
for the corresponding vertex. 

In Figure [3] cones in the fan of state polytopes for designs L\ and L2 
are depicted. As in each case the tridimensional cones form a partition 
of the first orthant, the figures show a slice of the cones when intersected 
with the standard simplex. The diagram for design L\ (left panel) shows 
27 cells, one for each model. The central symmetry of the diagram cor- 
responds to symmetry of models under permutation of indeterminates. 
Ordering vectors taken from the same cell will yield the same vertex (and 
corresponding model). Now in contrast with generic design L\, design 
1/2 produced the right panel in Figure [4] The diagram shows still some 
symmetry, but not central symmetry. This symmetry reflects the range of 
models computed for L2 in Example 1101 where only 11 models are iden- 
tifiable by L2, and ten models are in terms of X2 and X3. The following 
example illustrates changes in the fan by addition of one point to the 
design. 
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Example 15 Response surface design, central composite design. Con- 
sider the central composite design design in three factors built with axial 
points at distance y/2 and a full factorial design with points at levels 
±1. If no point is added to the origin, this design has 14 points and a 
combinatorial algebraic fan with 6 models. The models in the fan be- 
long to only two classes, one with monomials 1, z, z 2 , z 3 , z 4 , y, yz, yz 2 , y 2 , 
x,xz, xz 2 ,xy,xyz and the other class replaces z 4 by x 2 above. Addition 
of the origin to the previous design has a simplification effect in the fan, 
reducing to only 3 models, while it remains combinatorial. The only class 
of models is created by the list above together with the monomial x 2 . See 
depictions of both fans in Figure with the left panel depicting design 
without origin and the right panel after adding the origin. 




Figure 5: Algebraic fan of central composite designs in Example 1151 



7 Other topics and references 

The algebraic method in the form discussed here can be said to have 
started started with [44], |14j . [16] and the basic ideas were presented 
in the monograph [41] . A short review is [47] • More extensive work on 
the computation of universal Grobner basis with zonotopes appears in 
PQ. Applications to designs appear in mixtures [21] and [32] . Industrial 
applications were performed, perhaps surprisingly early: [27] [42| . 

For an excellent summary of the wider work in the field of Algebraic 
statistics: see [17] ■ One important topic omitted from this chapter, but 
important for conducting exact conditional test for contingency tables 
via Markov Chain Monte Carlo is the construction of Markov bases; see 
[26] [25] [24] [5]. Important applications to biology, which continue, are 
covered in [35] . Related and of considerable recent interest is the algebraic 
study of boundary Exponential models: [35] [5] [18] . 

Recent work showed the link between minimal aberration models and 
the border description of models in terms of Betti numbers of monomial 
ideals [53], see also extensive references in that paper. 
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